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Using series expansion based on the flow equation method we study the zero temperature prop- 
erties of the spin— 1/2 planar pyrochlore antiferromagnet in the limit of strong diagonal coupling. 
Starting from the limit of decoupled crossed dimers we analyze the evolution of the ground state 
energy and the elementary triplet excitations in terms of two coupling constants describing the inter 
dimer exchange. In the limit of weakly coupled spin— 1/2 chains we find that the fully frustrated 
inter chain coupling is critical, forcing a dimer phase which adiabatically connects to the state of 
isolated dimers. This result is consistent with findings by O. Starykh, A. Furusaki and L. Balents 
(Phys. Rev. B 72, 094416 (2005)) which is inconsistent with a two-dimensional sliding Luttinger 
liquid phase at finite inter chain coupling. 

PACS numbers: 75.10.Jm, 75.50.Eg, 75.40.-s, 78.30.-j 
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INTRODUCTION 



Geometrically frustrated quantum magnets above one 
dimension (ID) have attracted strong interest in recent 
years Q. From an experimental point of view this is 
due in part to the discovery of materials like the two- 
dimensional (2D) spin-dimer system SrCu2 (603)2 Q 
and the three-dimensional (3D) tetrahedral compounds 
Cu2Te205X2 with X=C1, Br which exhibit strong 

frustration of the magnetic exchange. From a theo- 
retical point of view, the prime interest in such sys- 
tems stems from the possibility of realizing quantum 
disordered states and spin liquid behavior. Particularly 
promising candidates in this direction are systems which 
comprise tetrahedra coupled into 2D or 3D networks. In- 
deed, experimental observations on pyrochlore materials, 
like CsNiCrFg |E or Tb2Ti207 suggest the pos- 
sibility of liquid-like behavior. In contrast to the lat- 
ter systems, which are essentially classical due to Ising 
anisotropy and spins larger than 1/2, the SU{2) invari- 
ant spin-1/2 Heisenberg antiferromagnet on the 3D py- 
rochlore lattice is an open issue. Experiments on the 
5 = 1/2 system Y2lr207 seem to indicate the pres- 
ence of quantum spin liquid state. To reduce the com- 
plexity of the 3D system, a 2D projection, i.e. the planar 
pyrochlore, also known as checkerboard or crossed chain 
model (CCM), which is shown in Fig^for Jp = Jd = Jo, 
has been considered as a simplified starting point. The 
CCM model and some of its anisotropic generalizations 
have been studied by employing a variety of techniques. 
Analytic methods include semiclassical S ^ 1 analysis 
|1EEE3> large-N limit of the Sp(N) model [iip,|li, 
series expansion (SE) by unitary transformation |lfil [ig 
and quantum field theory (QFT) methods 0, Il8l|. Be- 
sides, this model has been studied extensively by means 
of exact diagonalization (ED) techniques 0, [20I l2l| . 
These investigations have given a fairly complete picture 
of the phases of the system at different sectors of param- 
eters space. At the checkerboard point ( Jp = = Jq) a 
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FIG. 1: The dimerized planar pyrochlore lattice. Spin— 1/2 
moments are located on solid circles. Crossed dimers at site 
x,y (bold solid lines) are locally coupled via Jp (dashed lines) 
and non-locally coupled to nearest neighbors dimer by Jp and 
Jd (doted lines). On dimer exchange Jo will be set to unity 
hereafter. 



valence bond cr ystal (V B C) g round state with a spin gap 
has been found |l5llla . l20ll2ll |. For = Jq ^ model of 
Fig^tends to the 2D spin— 1/2 Heisenberg model on the 
square lattice, which is known to have gapless magnon 
excitations and AFM-LRO. From this limit, by means 
of a phenomenological approach within the framework 
of the Landau-Ginzburg- Wilson (LGW) theory, it has 
been suggested recently that the transition between the 
(AFM-LRO) phase and the VBC, as Jd(= Jq) are in- 
creased, could happen within a coexistence region or by 
a first order transition [l^ but this point deserves more 
investigation. Previous results from ED up to 36 spins 
for Jd = Jo in ref.ji^l predict a finite gap at Jd/ Jp — 0.85 
and an extrapolated closure of the gap at Jd/Jp — 0.65. 
Furthermore, semiclassical approaches predict stability of 
the AFM-LRO phase, in the case of S" = 1/2 at = ^0 
for Jd/Jp < 0.76. 

In the opposite limit of weakly coupled spin— 1/2 
chains, i.e. Jp <^ Jd = Jo in Fig^ there are contro- 
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versial results. On the one hand, previous findings based 
on QFT methods 0| suggested that the model exhibits 
a 2D spin liquid ground state, being an example of the 
so called sliding Luttinger liquid (SLL) which is charac- 
terized by the absence of LRO and the elementary exci- 
tations being massless, deconfined spinous. It has been 
suggested that the prime reason for this behavior stems 
from the momentum dependence of the completely frus- 
trated inter chain coupling. The latter vanishes exactly 
at the AFM nesting vector, suppressing a conventional 
magnetic instability. These results were corroborated 
subsequently by numerical studies [l^, where a range 
Jp/Jd G [0,— 0.8] for the extension of an effective ID 
behavior was suggested. 

On the other hand, by employing bosonization and 
renormalization group methods, Starykh et al. |l§l | have 
recently suggested that instead of exhibiting a SLL phase, 
the model of Fig^ enters a gapped crossed dimer phase 
for Jp <^ Jd = Jq, with a spin gap A which opens ac- 
cording to A/ Jrf ~ 0.675J2. 

In this situation it is desirable to obtain results from 
complementary methods, which will be the primary goal 
of this work. Therefore, in the following analysis we will 
employ a SE in terms of Jp and Jd as in Fig^ and will 
evaluate the ground state energy and the elementary ex- 
citations. Our findings will support the existence of a 
crossed dimer phase, in agreement with the results of 
ref.[l3. 

To perform a SE, applicable in the vicinity of the an- 
ticipated SLL phase we decompose the planar pyrochlore 
lattice into a set of crossed, dimerized chains as shown 
in FigH I.e. the model represents a 2D array of AFM 
coupled crossed spin— 1/2 dimcrs. We will set Jq — 1 
hereafter. The Hamiltonian reads 



H = Ho + V; Ho 



Y,ho; V = Y,iVp + Vd), (1) 



in which Hq refers to the local crossed dimer unit at site 
x,y, given by 

h„ = (Si- S3 + 82-84) (2) 
+ Jp(Si • 82 + 82 • S3 + S3 • 84 + 84 • Si) 

J p 

y 



,^13 + ^24 ~ 3) 



o V'^1234 ^13 '^24) 1 



where all spins are at site x, y and Si...„ = 81 + . . . + S„. 
Vp and Vd in Eq.lQ) represent the couplings between the 
crossed dimcrs at site x,y with their next nearest neigh- 
bors via Jp and Jd respectively, i.e. 



Vp — Jp[^4;x,y ' ^lix + l.y + ^3ix,y ' 82;2;+l,y 
+ S2;x,y • 8i:j,^j,+ i + 83^2:, j, ' S4-^x,y+l]', 



(3) 



Vd — Jd[^3;x,y ' ^l;x + l,y+l + ^2;x,y ' ^4;x-l,y+l]- 

Table m depicts the eigenstatcs of Iiq. From Eq.lPJ it is 





51234 


Sri 


^24 


q = E + 3/2 


l<So) 














\%) 


1 


1 





1 




1 





1 


1 


\St) 





1 


1 






1 


1 


1 




\Qt) 


2 


1 


1 


2 + Jp 



TABLE I: Eigenstates of a local crossed dimer unit ho at site 
x,y. Each state is labelled by the quantum numbers: S1234, 
S13 and 5*24 and the energy E. 



clear that each state of ho can be labelled by the fol- 
lowing quantum numbers: total spin 6*1234, the spin of 
the dimer on the positive (negative) diagonal, i.e. the 
'S'i3(S'24) and the energy E. As one can observe, for 
< Jp < 1 the ground state is \So), i.e. both dimcrs 
are singlets. Besides, for < Jp < 1/2 the first excited 
states are |7^)(|7^)) i.e. one of the dimcrs on the posi- 
tive(negative) diagonal is a triplet while the one on the 
other diagonal is a singlet. The other states arc: a total 
spin singlet |5t), a triplet |7t) and a quintet \Qt), all of 
them consisting of triplets on both diagonals. Only the 
energies of the latter three states depend on Jp. 



SERIES EXPANSION BY CONTINUOUS 
UNITARY TRANSFORMATION 

In this section we briefly describe the SE expansion in 
terms of V. To this end we rewrite Eq.Q as 



H = Ho{Jp = 0) 



N 

JpOp,+ E (JpOp^ + JdOdJ . (4) 



-TV 



Ho{Jp = 0) has an equidistant ladder spectrum (Table 
^1 which we label with the total particle-number oper- 
ator: Q = J2xyl('^p ^ ^^^'^ vacuum corresponds 
to the (5 = sector: |0) = Yix y\^0;x,y)- One-particle 
states correspond to the Q = I sector, whose basis states 
\Tfj.-x ,y) if- = PjIT') are composed of a local triplet on the 
positive or negative diagonal excited from the vacuum: 

\T,,x',y') ^ \'r^;x',y')^Ux,y^x',y' \So;X,y). The Q > 2 

sector is of multi-particle nature. 

The Opg operators in Eq.Q) refer to the local con- 
tributions proportional to Jp in Hq. Finally, the third 
term in the same Eq. represents the inter dimer coupling 
via V. The Op„ and Od„ operators ^ non-locally cre- 
ate(destroy) n > 0(< 0) energy quanta within the ladder 
spectrum of iJo( Jp — 0) . Note that only terms contain- 
ing two dimers at different x,y sites are present in V. 
This imposes the constraint N <2. 

Models of type Eq. I^J allow for SE using the flow equa- 
tion method of Wegner There, by means of a 
continuous unitary transformation, one maps H —>■ H^f / 



3 



with 



Ho 



Ef i jr — s js 



(5) 



r,0<s<r 



in which the Cr.s operators involve r — s and s products 
of type-Op„ and Od^ operators, respectively. We refer to 
ji^ l for further details. 

The main advantage of the flow equation method is 
that unlike if, the effective Hamiltonian H^ff is con- 
structed to conserve the total number of particles Q of 
Ho{Jp = 0) at each order in the expansion. This allows 
for the SE of several quantities from the bare eigenstates 
of Ho{Jp = 0) as we show in the following sections. 



GROUND STATE ENERGY 

The ground state energy allows for an important con- 
sistency check of our calculation with respect to exist- 
ing literature and other SE results. As mentioned above, 
there is evidence for a quadrumerYBC phase of the CCM 
at the checkerboard point (Jp = Jd = 1). This is based 
on ED |2ll| and SE starting from the quadrumer limit of 
the CCM [ilEl- Moreover, Starykh et al. H^, have 
predicted a first order transition or a coexistence phase 
at some intermediate value of < Jp < 1 and Jd — 1 
between the VBC and a crossed dimer phase. Therefore, 
in order for the present SE to give results consistent with 
them for Jp and Jd in the vicinity of the checkerboard 
point, we expect the ground state energy of the crossed 
dimer state to be larger than that of the quadrumer VBC 
at Jp = Jd 1. 

For < Jp < 1 the ground state of Ho{Jp — 0) is the 
Q = sector, i.e. |0) = Yl^ y \So)x,y Q-conservation of 
Heff allows for the calculation of the evolution of this 
state as a SE( Jp, Jd) by means of 



\Heff\0). 



(6) 



The actual calculation can only be done on finite clus- 
ters. To obtain results to 0{Jp,JJl), vaUd in the ther- 
modynamic limit, i.e. for infinite sized systems, ma- 
trix elements given by Eq.@ have to be evaluated on 
clusters with periodic boundary conditions (PBC), suffi- 
ciently large to avoid a wrap-around at graph length n 
[isf . In our case, the size of the Hilbert space, considering 
both perturbations Jp and Jd, restricts the calculation at 
n = 5 due to memory limitations. We find this enough to 
quantitatively describe the combined effects of Jp and Jd 
in the range of interest. At 0{Jp,J^) the ground state 
energy per spin and PBC, eg{Jp, Jd), is given by 
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FIG. 2: Ground state energy eg{Jp,Jd) per spin vs Jd 
for different values of Jp, Eq.(|7|l. As observed eg{Jp,Jd) 
monotonously decreases as Jp is increased from to 1, for 
fixed values of Jd- Eq.(|7J at the point of decoupled uniform 
spin— 1/2 chains, i.e. eg{Jp = Q,Jd = 1), differs « 1% from 
the exact value 1/4 — in 2 from Bethe ansatz. At the checker- 
board point Eq.0 gives eg{Jp = 1, Jd = 1) ~ —0.48 which is 
above the value obtained by employing a SE starting from a 
different point in parameters space (see text). This is consis- 
tent with the stability of the quadrumer VBC in the vicinity 
of the checkerboard point which has Cg ~ —0.51. 
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(7) 



As a check, note that eg(Jp ~ 0, Jd) reduces to the 0{J^) 
dimer expansion of t he g round state energy of the dimcr- 
ized spin— 1/2 chain l24[ l26l . 

From Bethe ansatz [27| the ground state energy for the 
homogeneous spin— 1/2 chain is known: Jp = 0, Jd ^ 
1) = 1/4 — In 2 I23. Comparing it with eg(Jp = 0, J^ = 1) 
from Eq.O we obtain a relative error of ^ 1%, which 
gives a rough estimation of the convergence . 

FigUl depicts the evolution of the ground state energy 
as given by Eq.© vs Jd for different values of Jp. As 
shown, eg{Jp, Jd) monotonously decreases as Jp is incre- 
mented from to 1, for fixed values of Jd- The posi- 
tion Jd,„ of the maximum of eg{Jp, Jd) vs Jd shifts from 
>/d„ = to Jd,-^ ~ 0.8 for Jp taking fixed values from 
to 1, respectively. At the checkerboard point our SE 
yields eg(Jp = 1, J^ = 1) ~ —0.48. Note that our value 
is larger than the value of « —0.51 which has been 
obtained from the quadrumer SE in ref.0|. While the 
expansion parameters and their respective orders cannot 
be cornpared directly between the present SE and that 
of ref . [l5j these findings are at least consistent with the 
stability of the quadrumer VBC in the vicinity of the 
checkerboard point. This, together with the fact that 
Ho{Jp = 1) has a degenerated ground state (see Table 
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1), suggests that Eq.|(7I), at least at the checkerboard 
point, does not represent the ground state energy of the 
system anymore, which is Cg ~ —0.51. 

TRIPLET DISPERSION 

Next we evaluate the dispersion of the one-triplet 
states. To this end it is necessary to diagonalize H^jf 
in the Q = 1 sector of Ho{Jp), i.e. the subspace spanned 
by the vectors \T^-x,y),fJ- = {p,n}- First we switch to 
momentum space 

\Tf,:k^^ky) = ^ [exp [ixK + iyky) \T^,:,.y)] , (8) 

in which PBC is assumed and is the number of sites. 
Because of the degeneracy of the Q=l sector, i?e// is a 
2x2 matrix in the basis given by Eq.©, whose entries 
ET^i,v-k^,ky are 

ETfi,v;k^,ky = {Tf^-^k^,ky\Heff\Tiy-k^,ky) 

= ^ exp {ixk^ + iyky)] , (9) 

x,y 

where fi,v = {p, n} and 

C^.jy-x.y = {Tp,-^x,y\Hef f\Ty-Qfi) (10) 

are the hopping amplitudes from a i^— type of triplet on 
site 0,0 to a /i— type of triplet on site x,y. Evaluation 
of the c^.i/:a:.t/'s is facilitated by symmetry which leaves 
only a subset of them to be independent. To obtain the 
hopping amplitudes valid in the thermodynamic limit, 
the c^,i/;2:,y's arc calculated on finite clusters with open 
boundary conditions (OBC), and a cluster-topology suit- 
able to embed all the paths of length n connecting the 
site (0, 0) with (a;, y) at 0{n) of the perturbation. From 
this, the triplet dispersion is obtained as 

W(k, Jp, Jd) = Et{K Jp, Jd) - (0|He//|0), (11) 

where k = {kx,ky) and the second term refers to 
the ground state energy which has to be evaluated 
on the identical cluster as c^,,y;o,o- The triplet energy 
i?T(k, Jp, Jd) results from diagonalizing the 2x2 matrix 
of Eq.jni), leading to two excitation branches. 

We have calculated the triplet dispersion a;(k, Jp, Jd) 
up to 0{Jp,Jd) Upon rotating k = {kx,ky) 

ikp,kn) along the positive and negative diagonals, by 
means of kx^y = (fcp T ^n) /2, the two branches of Ea. Hll() 
have to reduce to the triplet dispersion of the dimerized 
spin— 1/2 chains along positive and negative diagonals 
for Jp = 0. For consistency, we have checked that this 
is satisfied up to the order of our SE by comparing with 
existing SE results on the single dimerized chain [24 . l26l| . 



TRIPLET GAP 

Using the results from the previous section wc may 
now turn to one of the main objectives of this work, 
i.e. to study the evolution of the triplet gap under the 
combined effect of Jp and Jd- This will enable us to 
conclude on the stability of the dimerization along the 
chains and therefore to discriminate between the contro- 
versial proposals on the stability of the SLL. The triplet 
gap: A(Jp, Jrf) = mina'(k, Jp, Jd) is obtained from the 
minimum in k-space of the two one-triplet branches. For 
Jp < 0.5 we find the triplet gap to be at kp — kn — 0. It 
is given by 

A( Jp, Jd) = i- \jd - + ^4 + § Jd Jp' - ^4 

384^ 384 32 P 128 ^ 12288 

1309 .3 .2 _ 3895 2 3 6253 4 _ 2065 5 

^1536 P 4608 P ^ 3072 '' P 4608 J'' ^ ^ 

We have analyzed A(Jp, Jd) both by its plain SE and 
by using a biased-Pade approximationji^. The latter 
consists of an extrapolation biased by the known analyt- 
ical form of the gap-closure of the dimerized spin— 1/2 
chain as a function of Jd HJl, i.e. A(Jp — 0,Jd 
1) oc (1 — Jd)'^^^- Using this, Jd is replaced in favor 
of Jd = 1 - a;3/2 in Eq.JEI). Then, standard Fade 
approximationjs^ is performed around a; = 1 and finally 
X = (1 — Jrf)^/^ is re-inserted, arriving at a transformed 
Fade approximation in terms of the original variable Jd- 

Fig|2|depicts the evolution of A( Jp, Jd) as a function of 
Jd for different values of Jp. Solid lines show the results 
from the plain series, i.e. Ea. H12|l and dashed lines exhibit 
the biased-Fade (2, 2) approximants. As the inset of FigO 
shows, the approximant describes the closure of the gap 
at Jp = 0, Jd = 1 very well. This justifies the use of this 
extrapolation technique and shows that our expansion 
can be employed to analyze the coupling of the spin— 1/2 
chains via Jp. 

As a first result from Fig|3| both the plain series and 
the biased-Pade clearly show an opening of the gap at 
Jd = 1 for increasing values of Jp. As it is also obvious 
from these plots, this tendency persists for all Jd > 0.2. 
Second, Fig|31 shows, that the biased Fade approximant 
leads to observable differences from the plain SE only in 
the vicinity of the gap-closure. 

Now, wc consider the opening of the spin gap at Jd = 1 
for Jp <C 1. To this end, we first shift the zero of energy 
in order to suppress the finite gap of the SE at Jp=0 and 
Jd=l, i.e. we consider A( Jp, Jd = 1) — A(Jp = 0, Jd = 1). 
By this we account for the inability of the SE to provide 
for a truly gapless spectrum in the limit of decoupled 
AFM Heisenberg chains. As it is obvious from Figl^l 
this shift is very small for the biased Fade, i.e. « 0.003 
and PS 0.081 for the plain SE. The resulting spin gap 
versus Jp is shown in Fig0|whcre the solid line represents 



5 




FIG. 3: Triplet gap A vs Jd for increasing values of Jp, Ea. (ll2ll 
(solid lines) and biased-Pade(2,2) (dashed lines). As the inset 
shows, the biased-Pade(2,2) describes the closure of the gap 
at {Jp = 0,Jd = 1) very well. The effect of Jp opening the 
gap is evident from the arrow. 



the plain series and the dashed line corresponds to the 
biased-Pade(2,2) approximant. 

These results may be compared with recent findings of 
Starykh et al. Using QFT methods, i.e. renormal- 

ization group and bosonization, they have obtained an 
analytic expression for the spin gap, which is quadratic 
in Jp, namely A{Jp ^ 1, = 1) = cJp where c « 0.675. 
This function is shown in the same Fig^ with a dash- 
dotted line. Obviously, there is a very good agreement 
up to Jp « 0.1. This is remarkable, given the completely 
different nature of our approach and that of ref.[lg. We 
note that the same agreement is found when we em- 
ploy standard Dlog-Pade approximants instead of biased- 
Pade. Therefore, we are strongly tempted to claim that 
the stable phase of the CCM at finite Jp and at J^ = 1 
is indeed gapped. This is inconsistent with a SLL state. 

While Fig01 seems to suggest an ever-increasing spin 
gap, one has to remember, that at Jp = 0.5 the nature of 
the first excited state of /lo(Jp) changes from the triplet 
\Tp_n) to the singlet \St) (Table 1). This imposes an upper 
limit on the range of validity of the SE. 

CONCLUSIONS 

In summary, using series expansion, we have studied 
the zero temperature properties of the planar pyrochlore 
quantum magnet focussing on its behavior close to the 
crossed chain limit. In particular, we have considered 
the ground state energy, the elementary triplet dispersion 
and the spin gap. 

In relation with the ground state energy calculations 




FIG. 4: Evolution of the triplet gap A(Jp, Jd = 1) - A( Jp = 
0, Jd = 1) vs Jp starting from the limit of decoupled spin— 1/2 
chains {Jp — 0,Jd ~ 1). Plain series. Eg. 11211 . (solid line) 
and biased-Pade(2,2) approximant (dashed line). The shift 
A(Jp = 0, Jd = 1) has no effect on the functional dependence 
on Jp and it is made for a better comparison with the QFT re- 
sults by Starykh et al. [l|l, i.e. A(Jp < 1, Jd = 1) ~ 0.675Jp , 
(dash-doted line). Very good agreement between both tech- 
niques, up to Jp « 0.1, is observable. The same agreement is 
found by employing standard Dlog-Pade technique. 



our findings are consistent with previous results regard- 
ing the stability of the quadrumer VBC at the checker- 
board point. 

Our main result regards the existence of crossed dimer 
phase. Our calculations univocally show an opening of 
the gap A(Jp <^ l,Jd = !)• Moreover, our expan- 
sion (Ea. H12|l 'l. allows to route the physical origin of the 
crossed dimer phase by showing that it is adiabatically 
connected with decoupled crossed dimer units. 

On the other hand, the very good agreement with a 
recent field theory approach by Starykh ct.al. even 
for Jp < 0.1 consolidates the hypothesis about the exis- 
tence of a crossed dimer phase in the region of parame- 
ters space Jd = 1, Jp ^ 1 and indicates the absence of 
Luttinger liquid behavior, the so called sliding Luttinger 
liquid phase, in this region of parameters space. 
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